a+(980)-resonance production in the reaction pp — ► dir + r) close to the KK threshold 
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The reaction pp — > d-K+n has been measured at a beam energy of T p —2.65 GeV (p p =3.46 GeV/c) 
using the ANKE spectrometer at COSY- Jiilich. The missing mass distribution of the detected dir + 
pairs exhibits a peak around the n mass on top of a strong background of multi-pion pp — > dn + (rnv) 
events. The differential cross section d 4 a / 'df2 £ jdfi 7r +dp t idp 7r + for the reaction pp — » dir + n has been 
determined model independently for two regions of phase space. Employing a dynamical model 
for the a + production allows one then to deduce a total cross section of a(pp — » da + d-K + n) = 
(1.1 ± 0.3 s tat ± 0.7 sys ) fib for the production of K + n via the scalar a + (980) resonance and a(pp — > 
d-K+rf) — (3.5 ± 0.3 s tat ± l-0 sys ) /J,b for the non-resonant production. Using the same model as for 
the interpretation of recent results from ANKE for the reaction pp — > dK+K°, the ratio of the total 
cross sections is a(pp -> d(K+ K°) L=Q ) / 'a(pp -> da + -> d-K+n) = 0.029 ± 0.008 sta t ± 0.009 sys , which 
is in agreement with branching ratios in the literature. 



I. INTRODUCTION 

One of the primary goals of hadronic physics is the 
understanding of the internal structure of mesons and 
baryons, and their production and decay, in terms of 
quarks and gluons. The non-perturbative character of 
the underlying theory — Quantum Chromo Dynamics 
(QCD) — hinders straightforward calculations. QCD 
can be treated explicitly in the low momentum-transfer 
regime using lattice techniques which are, however, 
not yet in the position to make quantitative statements 
about the light scalar mesons. Alternatively, QCD- 
inspired models, which employ effective degrees of free- 
dom, can be used. The constituent quark model is one 
of the most successful in this respect (see e.g. 2]). This 
approach treats the lightest scalar resonances ao//o(980) 
as conventional qq states. 

However, more states with quantum numbers J p =0 + 
have been identified experimentally than would fit into 
a single SU(3) scalar nonet: the /o(600) (or er), / (980), 
/o(1370), / (1500) and / (1710) with 1=0, the k(800) 
and X*(1430) (1=1/2), as well as the a (980) and 
a (1450) (1=1) 3]. Consequently, the a // (980) have 
also been associated with KK molecules Q or compact 
qq-qq states 5]. It has even been suggested that a com- 
plete nonet of four-quark states might exist with masses 
below 1.0 GeV/c 2 Q. 

The first clear observation of the isovector ao(980) res- 
onance was achieved in K~p interactions |7j , and in sub- 
sequent experiments it has also been seen in pp annihila- 
tions 0, in n~p collisions 0, and in 77 interactions |lfjj| . 
Experiments on radiative ^-decays [ill have been 
analysed in terms of the ao / /o production in the decay 



chain — > 7CI0//0 — ► 77r°7y/7r 7r°. In pp collisions the 
a (980) resonance has been measured at p p = 450 GeV/c 
via /i(1285) — > flo 7rT decays and in inclusive mea- 
surements of the — > dX + reaction at p p = 3.8, 4.5, 
and 6.3 GeV/c [l4(. Despite these many experimental 
results, the properties of the ao(980) resonance are still 
far from being established. The Particle Data Group 
gives a mass of m ao = (984.7 ± 1.2) MeV/c 2 and a width 
of T ao = (50 - 100) MeV/c 2 3]. The main decay chan- 
nels, ttt] and KK, are quoted as "dominant" and "seen" 
respectively. 

An experimental programme has been started at the 
Cooler Synchrotron COSY-Jiilich [13 aiming at exclu- 
sive data on the ao//o production from pp, pn, pd and 
dd interactions at energies close to the KK threshold ^(| . 
The final goal of these investigations is the extraction of 
the ao //o-mixing amplitude, a quantity which is believed 
to shed light on the nature of these resonances [ItI Il8| . 
As a first stage the reaction pp — * dK + K° has been 
measured at T p =2.65 GeV, corresponding to an excess 
energy of Q=46 MeV above the K + K° threshold, using 
the ANKE spectrometer 19] . The data, which have been 
decomposed into partial waves, show that more that 80% 
of the kaons are produced in a relative s-wave, corre- 
sponding to the Oq channel [2(| • In this paper we report 
on the analysis of the reaction pp — > dn + ri, which was 
measured in parallel. 

II. THE ANKE SPECTROMETER AND DATA 
ANALYSIS 

Fig- IH shows the layout of ANKE. It consists of three 
dipole magnets (Dl - D3), installed at an internal tar- 
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get position of COSY. Dl and D3 deflect the circulat- 
ing COSY beam from and back into the nominal or- 
bit. The C-shaped spectrometer dipole D2 separates 
forward-going reaction products from the COSY beam. 
The angular acceptance of ANKE covers |#h| < 10° 
horizontally and |$ v | < 3° vertically for the detected 
deuterons {p d > 1300 MeV/c), and < 12° and 

I I < 3.5° for the pions. An H2 cluster-jet target [2lj . 
placed between Dl and D2, has been used, providing 
areal densities of ^5xl0 14 cm~ 2 . The luminosity has 
been measured using pp elastic scattering, recorded si- 
multaneously with the dir + data. Protons in the angu- 
lar range $ — 5.5° — 9° have been selected, since the 
ANKE acceptance changes smoothly in this region and 
the elastic peak is easily distinguished from the back- 
ground in the momentum distribution. The average lu- 
minosity during the measurements has been determined 
as L = (2.7±0.1 s tat±0.7 sys t)xl0 31 s™ 1 cm~ 2 , correspond- 
ing to an integrated value of L lnt = 3.3 pb" 1 at a proton 
beam intensity of up to ~ 4 x 10 10 . 
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FIG. 1: Top view of the ANKE spectrometer and simulated 
tracks of pions and deuterons from pp — > dn + rj events. 

Two charged particles, 7r + and d, were detected in co- 
incidence. Their trajectories, simulated with the ANKE- 
GEANT program package |22|, are shown in Fig. ^ 
Positively charged pions in the momentum range p^ = 
(600 — 1100) MeV/c were identified in the side detection 
system (SD) 0|23|, consisting of one layer of 23 start 
time-of-flight (TOF) scintillation counters, two multi- 
wire proportional chambers (MWPCs) and one layer of 
6 counters for TOF-stop. Pions were selected by a time- 
of-flight technique (Fig. The fast deuterons with 
momenta p d = (1300 - 2800) MeV/c (and elastically 
scattered protons with p p w 3400 MeV/c) hit the for- 
ward detection system (FD) [TsLI^ . which includes three 
MWPCs and two layers of scintillation counters. Two- 
dimensional distributions At vs. momenta of forward 
particles have been used for the deuteron identification. 
At was calculated as the time difference between the de- 
tection of a pion in the TOF-stop counter and a fast 
forward-going particle in the first layer of the FD scintil- 



lators. Two distinct bands corresponding to protons and 
deuterons are seen in this distribution (Fig. $2Jp). The 
selection of deuterons was achieved by cutting along the 
deuteron band. 
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FIG. 2: a) Time-of-flight of particles detected in TOF start 
counter 18 and stop counter 4; b) time difference At be- 
tween fast forward-going particles in the FD and -k + mesons 
as a function of the momentum of the forward particle. The 
dashed lines indicate the criteria for deuteron identification. 

The tracking efficiency for pions in the SD MWPCs 
was calculated as the ratio of particles in the proper TOF 
range with and without requiring a reconstructed track. 
The efficiency depends on the SD stop-counter number 
(i.e. 7r + momentum) and varies between 53% and 76%. 
For the FD MWPCs the efficiency has been determined 
for each of the six sensitive planes (two per chamber). At 
least two vertical and horizontal planes were demanded 
for track reconstruction and for the calculation of the 
intersection point with the remaining plane. Each plane 
has been divided in 20x20 subcells, and the efficiency 
of each cell has been calculated from the presence and 
absence of a hit in the reconstructed intersection. The 
average FD MWPC track efficiency for deuterons is 73%. 
The efficiencies of the scintillators and TOF criteria are 
larger than 99% |2^|. The efficiency correction is done 
on an event-by-event basis. 

III. RESULTS FOR THE REACTION pp -> dyr+ry 

The missing mass distributions mm(pp, d) and 
mm{pp, c?7r + ) for the selected dn + pairs are presented in 
Fig- El I n t ne (PP) dit + ) missing mass distribution a clear 
peak is observed around m(7/)=547 MeV/c 2 with about 
6200 events. The peak sits on top of a smooth back- 
ground from multi-pion production pp — » dit + (nit) (n > 
2). After selecting the mass range (530 — 560) MeV/c 2 
around the 77 peak, the missing mass spectrum mm{pp, d) 
exhibits a shoulder at 980 MeV/c 2 (Fig. EK, dotted), 
where the peak from the a,} (980) resonance is expected. 

Table [I] presents the differential cross sections for pp — > 
dn + ri for two regions of phase space where the acceptance 
of ANKE is essentially 100% for this reaction. Six vari- 
ables in the lab. system have been chosen to describe 
these rectangular areas: the vertical (9 y ) and horizon- 
tal (6 X ) angles and momenta of the two detected par- 
ticles, deuteron and n + . These angles are defined as 
tan(6» y ) = py/p z , tan(0 x ) = p x /p z - 
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FIG. 3: Missing mass distributions (a) mm(pp,d), (b) 
mm(pp,dir + ) for the reaction pp — > dir + X. The dotted 
histogram in mm(pp, d) (scaled by factor 6) corresponds to 
the selected area around the r\ peak (530 — 560 MeV/c 2 ) in 
mm(pp,dir + ) (indicated by arrows). 



d 4 cr/ dfJddf^-i- dp d dp 7r+ 
pb/(sr GeV/c) 2 


var 
0y, deg. 


iables, lab. sysi 
Ox, deg. 


em 

p, GeV/c 


71 ± 6 stat ± 20 sys 


d 

(-3°, +3°) 
(-4°, +4°) 


euteron variabli 
(-3.5°, +3.5°) 
pion variables 
(-11°, -3°) 


3S 

(1.4,1.6) 
(0.65,0.95) 


30 ± 4atat ± 9 sys 


d 

(-3°, +3°) 
(-4°, +4°) 


euteron variabli 

(-6°, -2°) 
pion variables 
(-11°, -3°) 


3S 

(1.8,2.7) 
(0.65,0.95) 



TABLE I: Differential production cross sections for the re- 
action pp — > d-K + rj. In both regions of phase space and 
non-resonant n + r/ productions contribute. For the momen- 
tum range pd = (1.4— 1.6) GeV/c the non-resonant -K + rj pro- 
duction should be dominant, because this momentum range 
corresponds to low masses of the n + rj system, where the 
production is suppressed. 



Unfortunately, due to the limited phase-space cover- 
age of ANKE, a partial wave decomposition of the type 
performed in Ref. [2(j is not possible in this case. Thus, 
in order to interpret the data in terms of and non- 
resonant Tt + r) production, we need to employ a model. 
This investigation will be described in the next section. 



IV. INTERPRETATION OF THE RESULTS 

The dir + r] final state can be produced via the a q res- 
onance, pp — > dciQ — > dn + r] (resonant production), or 
through the direct reaction pp — > dir + r] (non-resonant 
production). The production mechanism for the aj" has 
been studied theoreticallyin Refs. [U El El EI El 
According to Refs. [H |2(j, the cross section for 
production at T=2.65 GeV is expected to be ~ 1/ib 



while, for the non-resonant tt + t] production, different 
predictions exist, ranging from (0.6 — 1.4) fib [3(| and 
(1.6 — 3.3)^ib [U up to one order of magnitude more [26| . 

The differential cross sections from Table[|]contain con- 
tributions from both the resonant and the non-resonant 
reactions. We have performed a model-dependent anal- 
ysis using the calculated momentum distributions of the 
produced deuterons (see Figs.0Ji,b). The calculations are 
based on models describing the resonant process within 
the Quark-Gluon Strings Model (QGSM) (Fig. Et) 
and the non-resonant (Fig Eb) production via N*- and 
A-resonance excitation |29l l3ll |33j . The momentum dis- 
tribution for the resonant reaction is narrower, because 
low (ir + ri) masses are suppressed for a q production. The 
momentum distributions within the ANKE acceptance 
are shown in Figs. |3fc,d. 
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FIG. 4: Simulated acceptance for the deuteron momentum for 
the resonant (a,c) and the non-resonant (b,d) ir + ri production. 
Upper row: the initial model distributions |29l 1 3 1 L |32J, l33f : 
lower row: momentum distributions within the ANKE accep- 
tance. 
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FIG. 5: Diagrams describing the resonant ir +l q production 
within the QGSM |32| (a ) and the non-resonant via N* and 
A excitations [29ll3ll.l33]| (b). 

Within the models the non-resonant background has a 
negligible nr] s-wave contribution and thus does not in- 
terfere with the resonant amplitude. As a consequence, 
in the range pd = (1.4 — 1.6) GeV/c, the contribution of 
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the non-resonant process dominates whereas the resonant 
part is negligibly small (see shaded areas in Figs.0Ji,d). 
Thus, the cross section of the non-resonant contribution 
in this momentum range can be defined by fitting the r\ 
peak in mm(pp,dir + ) (Fig. Oi.) and extracting the num- 
ber of dir + r) events. The missing mass spectra have been 
described by the sum of a Gaussian distribution and a 
4 th order polynomial. 
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FIG. 6: Missing mass distribution mm(pp,dir + ) for the mo- 
mentum ranges pd — (1.4 — 1.6) GeV/c (a) and pd = (1.6 — 
2.8) GeV/c (b). 

In the momentum range pa = (1.6 — 2.8) GeV/c, where 
both the resonant and the non-resonant reactions con- 
tribute, it is possible to calculate the number of events 
from the non-resonant production taking into account 
the different acceptances. Then the number of o,q events 
is the difference between the total number of events un- 
der the r] peak (Fig. and the calculated amount of 
non-resonant Tr + r) events. 

Including all corrections, total cross sections <t& = 
(1.1 ± 0.3 s tat ± 0.7 S ys) A*b for the ciq production and 
Cn.r. = (3.5±0.3 s tat±l-0 S ys) A*b for the non- resonant ir + r) 
production have been obtained. It is obvious, that these 
numbers could change if different assumptions were made 
on the partial wave decomposition of the (non-)resonant 
contributions. However, we here refrain from investigat- 
ing further the uncertainty induced by the model assump- 
tions. 

Figure [7\i presents the results of GEANT simulations 
of the ANKE acceptance for multipion background (pp — > 
c?7r + (n7r) (n > 2)), and for resonant and non-resonant 
production of the dir + r] final state. The number of ini- 
tial events for each reaction is proportional to the known 
cross sections for multipion production [34| and the val- 
ues of the total cross sections for the resonant and non- 
resonant 7r + ?7 production given above. The shape of the 
simulated missing mass distribution mm{pp 1 d) including 
all channels is in good agreement with the experimental 
data (Fig.Eb). 

Data for the second decay channel, — > K + K°, 
have been obtained at ANKE in the reaction pp — » 
da^ — » dK + K°, simultaneously with the 7r + ?]_data. The 
measured total cross section is a(pp — > dK + K°) = (38 ± 
2 sta t ± 14 S ys) rib and the contribution of the (K + K°) L=0 
channel is ~ 83% |2jj. Assuming that s-wave (K + K°) 
production proceeds fully via the a J (980) resonance, 
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FIG. 7: a) GEANT simulations for multi-pion background 
(pp — ► dn + (nir) (n > 2)), and for resonant and non-resonant 
production of the dTr + ri final state in the ANKE acceptance, 
b) comparison between the simulated missing mass distribu- 
tion mm(pp,d) and experimental data. 

in accordance with the predictions of the models dis- 
cussed above, the ratio of the total cross sections is 
R = a(a£ -> (K + K°) L=Q )/a(a+ -> n+r]) = 0.029 ± 
0.008 sta t ± 0.009 sys . 
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FIG. 8: Ratio of the total cross sections cr(a,Q — * 
K + K°)/a(ag — > ir + r]) as function of Q. The error bar shows 
the statistical error only. 

For measurements at an excess energy Q » T ao ~ 
(50 — 100) MeV, the ratio of the total cross sections 
R = cr(a+ -> K+K°)/a(a+ -> tt+ rf) should not depend 
upon Q. However, when Q is of the order of r ao , the con- 
tribution from the — * K + K° decay channel _decreases 
more strongly due to the proximity of the K + K° thresh- 
old and the consequently limited available phase space. 
The calculated ratio of the two cross sections as a func- 
tion of Q is presented in Fig. [S] The calculation is nor- 
malised to R((KK)/{irr})) = 0.23 ±0.05, which has been 
measured at Crystal Barrel in the reaction pp — > aQir at 
Q = 768 MeV 35] . The parameters of the Flatte distri- 
bution 36] m = 999 ±2 MeV/c 2 , g m) = (324 ±15) MeV 
and r = g%jc/9^ri = 1-03 ± 0.14 are also taken from 
Ref. [3^. The errors of the Flatte parameters define the 
range of possible R values (shaded area in Fig. |SJ. Our 
result is in agreement with the calculated value. 
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V. CONCLUSIONS 

To summarize, data on the reaction pp — > dn+rj at 
T p = 2.65 GeV are presented. Differential cross sections 
for limited regions of phase space, corresponding to for- 
ward emission of the d and n + in the laboratory system, 
have been extracted. Dynamical models for the reaction 
were employed to obtain the total cross sections for the 
resonant n + r] production via the aj(980) and for the 
non-resonant TT + rj channel. For a$ production the value 
for the total cross section pp — > da^ — > dn + ri is in agree- 
ment with theoretical predictions based on the value for 
the pp — > d(K + K°)l = q reaction recently measured at 
ANKE |2(| and the branching ratio BR(KK/itt]) from 
the literature. 

Our results, together with those from ANKE on the 
decay channel Oq — > K + K°, are the first evidence for 
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